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Abstract 

We study the one-photon scalar pair production QED process on the expanding de Sitter space- 
time. Using perturbation theory, we obtain the transition probability and study its properties as a 
function of the expansion parameter uj. On flat space the process is forbidden by energy-momentum 
conservation. It is expected that for a dynamical background there is an energy exchange correlate 
to the strength of the gravitational field. We use momentum space plots and compute the mean 
production angle to illustrate this. We show that the mean angle grows with u, but also find 
that in the flat limit the fall-off is unexpectedly slow. To investigate this further we obtain the 
probability around different angular configuration, at leading order in m/ui, and find that the u 
dependence at small angles is very weak. We comment on the possible astrophysical implications. 
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I. INTRODUCTION 


In the absence of a complete theory of quantum gravity, the theory of quantum fields 
on curved spacetimes remains our best tool for investigating nature at the most funda¬ 
mental level. One of the most important results that have arisen from this approach is the 
phenomenon of gravitational particle production. In essence what this entails is that in a 
dynamic spacetime, pairs of particles can be produced purely out of the vacuum. Starting 
with the original paper by Schrodinger p], through the seminal work by Parker OT, up till 
the numerous present day studies, a considerable body of work has accumulated regarding 
the subject. Furthermore, the problem of self-interacting (scalar) fields on curved back¬ 
grounds has been investigated quite extensively. The main interest is in the renormalization 
on different backgrounds I5H7J and one of the important questions addressed is wether the 
different spacetimes are stable under perturbations (see for example 0E] and references 
therein). An adjacent topic, but which has been only sparsely investigated in the literature, 
is the effect of the dynamical background on the mutual interaction of quantum fields and 
inversely the contribution of the mutual interaction to particle production. Considering 
that in a physical setting the gravitational production is always accompanied by interaction 
processes, the study of the mutual interaction is imperative for obtaining a complete picture. 

Birrcll, Davies and Ford ra are the first to analyze the effects of held interactions upon 
particle creation in FRW universes. They argue that the interaction may significantly affect 
and, for bounce model universes, even dominate the pure gravitational production. The 
study of QED processes was undertaken by Lotze in expanding FRW universes mm, 
Tsaregodtsev et.al. in a radiation dominated universe mm and CotaescUjCrucean et.al. 
in de Sitter space [T6HT8] . Furthermore there are a number of papers dealing with the 
quantum held theoretic treatment of radiating charges in different FRW universes [T!M?T] . 
at leading order in the WKB approximation or the perturbation theory. 

The studies involving de Sitter space (dS) are particularly relevant due to the cosmological 
significance of the spacetime. In this paper we present a detailed analysis of the transition 
probability for one-photon scalar particle pair production on the expanding patch of dS. 
We work in the in-out formalism, with the S-matrix approach. In order to focus on the 
contribution of the interaction to the transition probability, we consider the hnal vacuum 
state to be identical to the in vacuum, and define the notion of particle with respect to the 
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in modes. The resulting amplitudes are sometimes called in-in amplitudes in the literature. 
The notion of particle is defined with respect to the in-modes. The resulting probability is 
interpreted also in terms of added-up probabilities, introduced in ref. [223 an< 4 described in 
Appendix [B} 

The paper is structured as follows. Section I contains the necessary ingredients for setting 
up the perturbation theory for the scalar QED on de Sitter space. In section II we obtain 
the transition probability and perform an extensive analysis as a function of the expansion 
parameter. In section III we summarize the main results. We work in natural units where 
h — c — 1 . 


II. PRELIMINARIES 


The expanding Poincare patch of the de Sitter spacetime is described by the following 
metric, written in both planar and conformal coordinates [23] : 

ds 2 = dt 2 — e 2uJt dx 2 (1) 

= (a :t c Y 2 (dt 2 c — dx 2 ), 

where u is the Hubble parameter and the conformal time is defined as ut c = —e _a; / t c G 
(-oo,0). 

We consider a massive scalar field minimally coupled to gravity. The Klein-Gordon (KG) 
equation in this geometry is: 

(< d 2 — e~ 2ult A + 3 ud t + m 2 ) (p(x, t) = 0 (2) 


The positive frequency plane-wave solutions of the KG equation are 


fA x ) = o 


g*W 2 JJ{1) ( P_ e -wt\ ipx 

\M 2tt) 3 / 2 v \u ) ’ 


( 3 ) 


where H^, are the Hankel functions with index v = p? — |, and m = up is the 
mass of the scalar field. These modes define the Bunch-Davies vacuum [2Blj . 

The Maxwell field is quantified in the Coulomb gauge 123. given by A 0 = 0, (y/^gA l ). t = 
0. Because the field is conformal, the mode functions can be straightforwardly obtained 
from the flat space solutions I23 : 


W 


X = 


—ikt c +ik x 


k ’ x (2vr) 3 / 2 a/2 k 


£\(k)e 


—2cot 


( 4 ) 
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where £\{k ) are the polarization vectors. 

With the above expressions ([3]) and Q, the usual mode expansion for the free scalar and 
Maxwell fields can be written as: 

<p(z) = j d 3 p (a(p)fp(x) + b\p)f* p {x )) (5) 

A \ x ) = J] [ d 3 k (c\(k ) w\ x (x) + c* x (k ) wgjz)') (6) 

A d 

A°(x) = 0, i = 1,2,3, 

where a\W and c t are creation operators for scalar particles, antiparticles and photons, 
respectively. 

We use a general prescription for interacting field theories [26] and consider QED processes 
in the first order of perturbations. The 1st order term in the perturbation expansion of the 
S matrix is given by: 

S (1) = -e f d A xyf^g(<p + (x) d M p{;x))A^{;x), (7) 

where \f—g = e 3ut and f d g = f(dg ) — g(df ) is the bilateral derivative. We have omitted 
the four-point interaction term, which is second order in the coupling constant, because it 
does not contribute to the processes studied here. 

The transition amplitude for these first order processes is defined as: 

Ai->f = (out |S' < - 1, | in). (8) 

With this, we have all the preliminaries in place in order to study tree-level processes. In 
our previous work [T8] we have considered the emission of photons by a scalar particle, and 
here we investigate the process of one-photon pair production. 


III. ONE-PHOTON PAIR PRODUCTION 

Within the framework introduced in the previous section, we consider the process of 
scalar particle pair creation by a single photon, illustrated by the Feynman diagram ([Tj) . 

In flat spacetime this process is forbidden at tree-level. This can most easily be seen 
by changing to the center of mass frame of the resulting pair of particles. Due to energy- 
momentum conservation, the momentum of the photon would then have to be vanishing 


4 



FIG. 1: Feynman diagram - pair production 

while at the same compensating for the rest energy of the massive particles, at the least. 
An analogous situation is instead the second order process of one-photon pair creation in an 
external field, commonly in the field of a nucleus [28] or in a strong magnetic field [29]. The 
first is relevant because of the possibility of experimental confirmation, while the second is 
relevant in astrophysical settings. In our case the background plays the role of the external 
(gravitational) field. 

This process is also of interest from another point of view. The interpretation of tran¬ 
sition amplitudes and more generally of what represents a measurable quantity on a non- 
asymptotically flat spacetime is still a matter of some debate. It is thus desirable to deal 
with quantities whose interpretation is insensitive to the definition of the out state for ex¬ 
ample. The concept of added-up probability, introduced by Audretsch and Spangehl |22j, is 
such a quantity. 

In this approach one calculates the probability for transition between two states of the 
Maxwell field with different occupation numbers, irrespective of what happens to the massive 
field. This is accomplished by summing over the complete Fock space of the scalar held in 
the final state. This method can be employed because the Maxwell held is conformal and 
thus the vacuum state is the same for all times, i.e. photons can not be created freely from 
the background. A brief summary of the method of added-up probabilities can be found in 
Appendix [B} 

For the specihc case of pair production by a photon, the added-up probability means the 
transition from the in-state containing a photon (an no scalar particles) to a state with no 
photons (vacuum of the Maxwell held), summed over all possible outcomes for the scalar 
held. In other words, the probability that a photon ’’disappears”, regardless of ’’into what”. 
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What is remarkable in this case is that the diagram (|Tj) is the only one that contributes 
to the added-up probability [12, [T3}. This is not true for the time-reversed process (see 
Appendix [B]) . The added-up probability can thus be used as a Rosetta stone to translate 
between in-in and in-out probabilities. 


A. Transition probability: expression 

The transition amplitude which describes the one-photon scalar pair production has the 
following expression: 

A(p,p', k) = (1 m , 1 ( p -v) | S w | 1 {%X) ) 

= -e / d 4 x^g (f£(x) di /^(z)) w^ x (x). (9) 

Introducing the solutions (|3]) for the scalar held and Q for the Maxwell held, the transition 
amplitude becomes: 

ie7r(p' — p) ■ e [k) 


A(p,pk) = 6 3 (p + p — k ) 


(27t) 3 / 2 v / 32/c 




( 10 ) 


where 


I^{p,p‘,k)= / dppHl a \ V 'ri)Hy<j> p)e ik \ 


( 11 ) 


The form ( [TT| for the temporal integral was arrived at by changing the integration variable 
to ?? = —t c . 

The tree-level (scalar) QED amplitudes on dS space all have the same structure, with a 
momentum-conserving delta function arising from the spatial integral as a consequence of 
spatial translation invariance, and with the temporal integral contributing the non-trivial 
physics. 


(2 2 ) 

The calculation of the integral iy is given in Appendix ( A 1). From the form (A7) of the 


amplitude we observe that it depends only on the ratio p = m/cu and the involved momenta. 
Unless stated otherwise, in the following we shall consider unit mass in the formulas and 
plots. 

We add to the integral ( p~T| an exponential factor e~ ev that acts as a switch-off for the 
interaction for large times, the decoupling time being of order 1/e. The e is a small positive 
parameter that reinforces the convergence of the integral. Effectively we consider it to be a 
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small but finite constant in our calculations, and discuss the vanishing limit at the end of 
the paper. 

The partial probability, averaged over the photon polarizations, is obtained as: 


'P(p,p',k) = ^ A(p,p’,k) 

Z A=±l 


( 12 ) 


The delta term can be handled, following the prescription from flat space scattering theory, 
by writing 

( 27 ,'fs 3 QTp )| 2 = (27r)«<s‘ 3 >(0)^ Qr?) = (2jr) 3 V 3 S :i (j^f) , (13) 

where V is the comoving volume. The physical volume is V p h ys = V e 3u}t = ( ujr]) 3 V . One 
then usually considers the probability per unit volume. In the following we will shorthand 
this quantity with just the probability, remembering that what we really mean is in fact the 
probability per unit comoving volume. 

The polarization term can be easily obtained by making use of the relations p' = k — p 
(momentum conservation) and k ■ £\{k) = 0. 


(p'-p)-£*$) 


2 p-£\(k) 

(p -k) 2 


(14) 


=4 r - 


k 2 


= 4 p 2 sin 2 9 
4 p 2 p' 2 sin 2 x 
= k 2 


where in (14) we recognize the projection of the momentum onto the plane defined by the 
polarization vectors (and perpendicular to the direction of k ), and 9 represents the angle 
between p, k, while x represents the angle between p,p'. 


Gathering all terms and using the notations from Appendix (A 1), the probability becomes: 

/t-,/ ->/ A r 3 /-* , ->/ e 2 7r 2 p ,2 p 2 sin 2 x f . 

V{p,p ,k) = &\p + p - k) —— (15) 


16 (27r) 6 fc 3 

x 1 9 v(p,rf,k) + g-„(p,p',k) + h v (p,p', k) + h- v (p,p’, k )\ 2 . 


The quantity (15) is to be interpreted as the momentum space distribution of the total 
probability: 

dP 


V{p,p,k) = 


d 3 p d 3 p' ’ 


(16) 
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from where the total probability is defined as: 

P(k) — I d 3 pd 3 p' V(p,p',k) 
We shall also consider the following the quantity 

V(k,p,6)= [ d 3 p' V(p,p',k). 


( 17 ) 


(18) 


Notice that the delta function makes the above integral trivial. 

For later use, we obtain the expression of the probability for the particular value p = y/2, 


where the Hankel functions have simple analytical forms (Al) and the temporal integral can 
be readily evaluated: 


I±i\p,p':k) = I dr)r) 


?-wv 


e ~iPV e ik V 


2 1 


7T a/ pp 7 
2 1 


np’p v Trp V 

dr]V e -i( P '+P-k)v 


(19) 


7T a /pp 7 p' + p — k 


The probability thus becomes: 

V(k,p,6) = 


9 

e p 


sin 2 6 


( 20 ) 


(27t) 6 4 k p' (p + p' — k ) 2 + e 2 
where we have restored the e parameter for clarity, and it is understood that p' = 
\Jp 2 + k 2 — 2pkcos9. The simple form of the probability in this case makes it particu¬ 
larly suitable for numerical and analytical manipulations, and will be exploited accordingly 
in the following. 

We note that for this particular case the mode functions reduce to the flat-space plane- 
waves, and the scalar held effectively behaves as a massless conformally coupled held. One 
would then expect that the probability be null as in the hat space case. The probability is 
indeed equal to the hat space one, but neither of them are in fact null. This can be seen by 


identifying the conformal time q with the Minkowski time in eq. 19 We observe that the 
integration is only over the semi-infinite interval. This is because the expanding patch of 
dS is only locally conformal to Minkowski space. The integration up to a finite time in hat 
space can be understood as a sudden decoupling of the helds at a finite time. This would 
then lead to transient effects giving the probability a non-null value. Translating back to 
the expanding dS space, what in hat space was due to (unphysical) transient effects, now 
represents the physical probability rendered so by the nature of the spacetime. 













B. Transition probability: analysis 


The probability can be evaluated analytically for m uj (weak gravitational regime), 


with the aid of an approximation. From the expression (15), the h± v terms which contain 
an Appell F 4 function are the ones that resist analytical manipulation. We have found that 
the approximation 


F 4 (l, 3/2,1 + v, 1 - v, x, y) ~ 1, 


( 21 ) 


works very well for large p. A more detailed discussion can be found in Appendix (A3). 


An alternative option is to numerically integrate the temporal integral in (10), which we 
shall do for the domain m ~ oo (strong gravitational regime), where the above approximation 


is not accurate. Finally, we will make use of the analytically evaluated probability (18) for 
the case p = \/2. 


The probability is best visualized in cartesian momentum space surface plots, with p x 
denoting the direction parallel to the photon momentum and p y the perpendicular. As can 
be seen from fig. ([2]), the probability takes significant values only in the interval p x € (0, k ) 
and at small angles (small p y component), falling abruptly for negative and large momenta. 
The overall conclusion to be drawn is that even tough there is the possibility of energy 
exchange with the dynamic background, the probability is peaked in the vicinity of the 
energy conserving case (9 — X —> 0, k ~ p + p'). To get an intuitive picture one can think of 
the configuration of momenta in Compton scattering in flat space. In that case the scattering 
at backward angles is forbidden by the simultaneous energy-momentum conservation law, as 
is the case where the photon gains momentum as a result of the process. In the somewhat 
analogous situation for our process of pair production in de Sitter space, the backward 
’’scattering” and production of large momentum pairs, although not forbidden, are heavily 
suppressed. On the other hand, a significant relative increase in the probability of production 
at intermediate angles, corresponding to production in the forward quadrants, is clearly 
visible as the strength of the gravitational field increases. Also, notice that while in a weak 
field the probability has a maximum at p ~ p' ~ k/2, in the strong field case the probability 
is the same throughout the interval (0, k ). 
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(a) n = V2 


(b) n = 2 
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(c) fi = 5 (d) /x = 10 

FIG. 2: Momentum space surface plots of the transition probability, in cartesian coordinates, for 
k=l and e = 10~ 2 . The x-axis is taken in the direction of the photon momentum k. The 
probability being invariant under rotations around the direction of k, the problem reduces to 

scattering in a plane. 


C. Mean production angle 


To characterize the angular behaviour of the probability with varying strength of the 
gravitational field, we calculate the mean production angle for various values of the expansion 
parameter. If we consider the transition probability as the distribution function for the 
variable 9 , we can obtain the mean production angle as: 


ff 6V(k,p, 9)p 2 dp sin9d9 
ff V(k, p, 9) p 2 dp sin OdO 


( 22 ) 
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The straightforward calculation of the mean angle is not possible however because the prob¬ 
ability has an UV divergence. This is highly unexpected, especially considering that in the 
flat limit the mode functions reduce to the familiar plane-waves, as was shown in ref. [33], 
and the probability vanishes accordingly. To trace the origin of this divergence we note that 
in the ultra-relativistic case (p ;§> m,prj 1), the mode functions reduce to the flat space 
plane waves, but with time coordinate p. One would expect that the amplitudes then re- 




p p 

(a) (b) 


FIG. 3: The momentum distribution of the probability for p = \/2 (a) and p = 10 (b), for k = 1 
and various angles. Note the sudden drop of the probability for p > k, and also the many orders 
of magnitude difference in the large p behaviour for the two cases. 


produce the flat space ultra-relativistic amplitudes for the analogous process in an external 
field, for example. As we have argued above however, because the range of the conformal 
time is restricted to the semz-inhnite axis, the amplitude will actually be similar to that in 
flat space, but where the interaction has been suddenly decoupled at time t = 0. This would 
then lead to transient effects and an UV divergence when integrated over the momenta of 
the particles, as was reported also in ref.[31j. 

In order to extract a quantity which gives at least qualitative information about the 
angular behaviour, we effectively integrate the momentum up to a maximum value p m ax- 
Observing that the probability is highest in the interval (0, k), as can be seen from fig.(j2]) 
and Q, we consider it sufficient to take p ma x ~ k in order to illustrate the variation of the 
mean production angle. Increasing p ma x accentuates further the increase of the mean angle. 
Listed below is a table with the results obtained for the mean angle, for different values of 
Pmax and e: 
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h 

e = 0.01 

Pmax — 2 k 

e = 0.01 

Pmax — 10& 

e = 0.001 

Pmax = 2 k 

e = 0.001 

Pmax — 10 k 

10 

6.76 

6.76 

3.53 

3.53 

5 

9.40 

9.54 

5.05 

5.37 

2 

20.40 

23.90 

25.40 

31.25 

V 2 

45.11 

48.45 

35.66 

60.40 


TABLE I: (0)°- mean emission angle 


As one can see, even tough the values change for the different parameters, for all cases 
an increase in the strength of the gravitational field (decrease in /i) leads to an increase 
in the mean production angle. We note that the two intermediate values may not be very 
accurate because of the errors in the numerical integration, but they are in conformity with 
the overall trend. 


D. Flat space limit 


As we have mentioned, this process is forbidden in flat space because of incompatibility 
with energy-momentum conservation. This implies that the probability in dS space must 
vanish in the limit to —)• 0. To get a better picture of the behaviour, we keep the leading term 
in p = m/u in the probability, and seek the approximate expression around three different 
angular configurations between the involved momenta. 


The main contribution to the probability in the flat limit comes from the gk term in (15) 


This is not easy to see however, because the /i-dependence is intertwined with the angular 
dependence. Using the approximate relations 


v ~ ip, sinh(7r/r) ~ cosh(7iyi) ~ - e 71 ^, 


the leading contribution to the probability becomes: 


16 (27t ) 6 k 3 

e 2 7r 2 sin 2 (x) 

4 (27r) 6 kpp' 


(23) 


(24) 


^e n,J ' 2 F\ ( - - f/i, - + f/i, 2, 


1 + cos(x) 




x 




1 - cos(x) 
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where we have denoted ((e) = ^ 7 . 

The remaining terms fall off as higher powers of e _7r/ h We have kept track of the e parameter 
because it plays a crucial role in the form of the probability, as we shall see below. 


Using the formulas from Appendix (A2) to approximate the hypergeometric functions 


contained in (24), we obtain the expression of the probability in the particular cases x — 
0, 7t/ 2, 7r. 


V(x — o) — 
V(x = 7T/2) ~ 
V(x — 7r) ~ 


X 


(27r) 6 kpp' x 4 + 16C(e) 2 


7T/i 


(27t) 6 kpp 1 2 


7 Vfl 


7r 






(25) 

(26) 
(27) 


(2ir) 6 kpp' (jr - x) 4 + 16C(e) 2 
The difference among the three particular angular configurations is striking. For the one- 
photon pair production process in dS space we have found that the probability falls of as 
e -«(x)/j., w fth a(x) > 0. This behaviour carries over to other tree-level processes, and most 
likely also to higher order processes. Most importantly, a —> 0 for x — ® 0, which means 

that small angle pair production is dominant in the small expansion parameter limit. 


By extrapolating (25) to a small vicinity around 9 = 0, we conclude that because here the 
/i dependence is weak, the probability can remain large even for a weak gravitational held. 
The probability for the process, being a function of the ratio /j = m/u, given a small enough 
mass for the scalar held, can be significant even for the present day expansion. Our results 
then suggest that this can still be the case even for particles with Compton wavelengths 
much smaller than the Hubble radius (m co). This could have potentially interesting 
astrophysical consequences. 


IV. CONCLUSIONS 

We have studied the hrst order QED process of one-photon scalar pair production, on 
the expanding Poincare patch of the de Sitter spacetime. 

The transition probability was evaluated with three different methods for three different 
domains of the expansion parameter: a) by approximating the analytical expression in the 
weak field case (m oj), b) by numerical integration for the strong gravitational domain 
(m ~ oo) and c) by analytical integration, facilitated by the simple form of the mode functions 
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for a particular strong field case (m = y/2u ;). The coherence of the results acts as check for 
the validity of the methods. 

The resulting probability distribution was represented as a surface in the momentum 
space (hg|2]). We have found that there is a moderate interaction (i.e. energy exchange) 
with the background: 

• in general the probability is concentrated around the configuration that is closest to 
the energy conservation condition (k ~ p + p' , x — 0 —> 0) 

• by turning up the strength of the gravitational held, the production at intermediate 
angles is enhanced, while remaining vanishingly small for backward scattered and high 
momentum pairs. 

To make the above description more precise, we have calculated the mean emission angle 
for various configurations and have found accordingly that it increases/decreases with the 
increase/decrease of the expansion parameter. In the hat limit the mean angle goes to zero 
and probability vanishes as expected. Surprisingly, the fall of the probability is rather mild 
for small angle production. 

To investigate this behaviour, we have approximated the expression of the probability 
in the vicinity of different angular configurations, while keeping only the leading term in 
p = m/uj. We have found that the the probability falls oh as with a(x) > 0. 

This behaviour carries over to other tree-level processes, and probably also to higher order 
processes. Most importantly, a — > 0 for x — @ ~* 0, which means that small angle pair 
production is dominant in the small expansion parameter limit. Noting that the probability 
is a function only of the ratio m/uo, given a small enough mass the probability can be large 
even for the present day expansion. Our results then seem to suggest that even for a scalar 
held with Compton wavelength much smaller than the Hubble radius (m ~ oj in natural 
units), the probability of pair production around small angles can remain non-negligible 
for the present expansion. The fact that we do not observe this process for example for 
the CMB photons, gives a lower bound for the possible masses of such ultra-light particles. 
Such scalar helds, with Compton wavelengths between Hubble and galactic order, have been 
popular in recent years in beyond the standard model physics, as candidates for dark matter 
and dark energy. 
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We signal that these small angle results have to be taken with a dose of skepticism. 
This is because the quantitative description is highly dependent on the switch-off parameter 
e. When the vanishing limit is taken the probability is divergent. This is not however 
a dangerous divergence, being rather a universal trait of theories where massless particles 
are involved. It arises also in flat space, for example, in the case of a radiating charge 
kept at constant acceleration by an external source or evolving in the held of a nucleus 
(bremsstrahlung). 

It would be interesting to see whether this divergence can be dealt with in a similar 
manner as in the hat space counterparts. Indeed the entire divergence structure of these 
transition probabilities is interesting and might have significant consequences. 

a: for the bremsstrahlung case, the (IR) divergence is eliminated when the higher order 
contributions and also the elastic scattering are summed. It would be interesting to 
see whether the e —> 0 divergence in the de Sitter amplitudes might also be eliminated 
when the higher order interactions and also the elastic scattering oh the gravitational 
held are taken into account. 


b: the (UV) divergence noted in section III C is more bothersome. It arises because of the 
finite (conformal) time integration in the amplitude, and we have argued that it is 
similar to the transient effects arising when one performs a sudden decoupling of the 
helds in a process on hat space. This might represent a breakdown of the applicability 
of the eternal dS space in physical contexts. 


c: the dS QED probabilities have an additional divergence in the m/u ; = {0,1.5} domain, 
stemming from r) —> 0 end of the temporal integral. This is the famous (IR) divergence 
of dS space which is relevant for inflation, and is still an open problem to date (for 
a comprehensive list of references see 0 )- It might be worth investigating if the 
interaction could have a word to say regarding this issue. 
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Appendix A: Mathematical Toolbox 

1. Bessel functions 


Here we will list the relations and integrals involving Bessel and Hankel functions, that 
enter in the calculation of the amplitude. [321 [33] 

For the particular value of the index ±| the expression of the Hankel functions is: 


2 

r<1 !( 

2 



2 

7 TZ 



r ( 2) 

2 

r (2) 
~ 2 



(Al) 


The Hankel functions can be expressed in terms of Bessel J functions as follows: 

J- V {z ) - e~ mv J v [z) 




l Sill 7 TV 


(A2) 


— e dyjz) 


l Sill TTV 

The temporal integral can be thus turned into 4 integrals of the form: 


J dq rj J± v {prj)J± v {p'ri) e lkq J dpr) J± v {prj)J^ v {p'rj) e lkr) . (A3) 

o o 

The first type of integral, which contains Bessel functions with equal sign indices, can be 
solved by using: 

OO 

/a rj 

dp V J± u (prj) J± u (p'rj) e ikr, - eri = ~ ^ ppl)3/2 ( z ) ( A 4) 

o 


d 

dz 


Quiz) 


7 jv{v + 1 ) 

4 sin nv 


2jFi 


^1 — Uj 2 + v] 2 ; 


+ 2 Ai 


^1 — v, 2 + v] 2; 


i+*y 

2 )\' 


(AS) 
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where z = p +v 2p p +l ^ ? an d the ± branch is selected according to the sign of -Im(z). The 
second type of integral, containing Bessel functions with opposite sign indices, is given by: 

OO 

J dr] 1 7 J v {jprj) J- V {jp'r)) e lkv ~ €TI = (A 6 ) 

o 


p\ u ( 1 \ 2 sin(7rz/) / 3 p 2 p' 2 

]y) \ik) 7 TV \ 1, 2 ,1 + U, ~ V ' (k + ie ) 2 ’ (k + ie ) 2 


With the above formulas, the temporal integral can be written as: 


e tnv lj?’ 2) = [gu{p,p', k) + g-u(p,p’, k) + h„(p,p ', k) + h-„(p,p', k)]., (A7) 


where we have introduced the notations 


g ±v(p,p\k) = 


=PZ7 TIP 


ik (v 2 ~\)e 
A(pp') 3 / 2 cosh(7ri/) sinh 2 (7r/c) 


ie* tm ' 2 F 1 I ± n, ^ =F i/; 2; 


3 

2 


h ± ^(p,p',fc) = 


k 


-2 


7rz/sinh(7rz/) 


F , 3 , 3 _ 1 + * 

2*1 ( 7^ ± G 2 T 2; ~T” 

p \ ±v f 3 p 2 p' 2 


(A 8 ) 


2. 2 A] hypergeometric function 

In order to find the asymptotic form of the probability, in the flat limit, we need to 
find an approximate form for the Gauss hypergeometric functions that are contained in the 
leading terms. More explicitly we wish to find the form of 2 Fi (§ — v, § + g2,r), for the 
distinct cases x = { 0 , |, 1 }, in the limit p —» oo, v —» ip. 

With the use of the relations [32] 

2*1 (a, P, 7, x) = (1 - x) 7_Q_/3 2 * 1(7 - a, 7 - P, 7,a) ( A9 ) 

2 Fi(a,ft 7 ,l)= Re(l)>Re(a + l 3), (A10) 

r(7-a) r (7-p) 


17 













we find: 


2 F 1 ( ^ - U ^ + ", 2, x ~ 1 j ~ (1 - x) 1 2 -Fi Q + ¥, ^ - ©, 2,1 


(AH) 


= (1 -x) 


-1 


r( 2 ) r(i) 


r (| + f/i) r (| - z/i) 
(1 — x) 1 


1 + z/i) r (i + z/i) (i - z/i) r (i - in) 

(1 — © 1 


(i + ©) | r Q + i/i) | 

1 COsh( 7 T/i) 


(!-®) ttQ 


For the case x — 1 we require the formula: 


2 -Fx ( 2cc, 2/3, a + /3 + -; — 2 ^^ 


^ — A 2 F[ ^a,/3, -]y^j + B yd/ 2 Fi fa + -,/3 + -, 


A = 


r (a + jd + 1) v/7r 


_ o* _ 3 , ¥ 

P ~4 + Y- 


r (“ +1) r (^ +1) ’ 

In our case y = 0, and the right-hand side of ( |A12 ) reduced to A. 
Using also the limit: 

lim |r(a + i/3)| e~^ Z?© - ") = V^n, 


(A12) 


(A13) 


we obtain: 


2 F 1 ( ^ - v, ^ + u, 2, x = /*- 1 = 


7 T 


r (! + f) r (!-f) 


7 T 


Q-f)r(! + f)(i + f)r(!-f) 


7 r 


l(l+/i 2 ) 27 r(f) 1/2 e T 



4 / 2 1 


Finally: 


^uu+if© 1 '' 2 - 


2 F 1 I | - V, / + v, 2,0 I =1 


(A14) 


(A15) 
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3. Appell F 4 function 


From the definition of the Appell F 4 hypergeometric function, we have: 


OO 

F 4 (l,3/2,l + u,l-u,x,y) = ^ 

m,n =0 


0) m+n (3/2) m+n x m y n 
(1 + u) m (l - u) n m\ n\ ’ 


(A16) 


where the Pochhammer symbol means (a) m = a(a + l)...(a + m — 1) = T(a + m)/T(a). 

We are interested in the limit fi — > oo, v — > i/i. As the parameter fi appears only in 
the denominator, it is reasonable to consider the approximation F 4 ~ 1. This is by no 
means obvious, because our function does not satisfy the absolute convergence criterion: 
Vx + y/y<l. 

A comparison of the approximated against the numerically evaluated probability can be 
seen in fig [4j We have found that for m > 5co this approximation is very good. However, 
increasing the switch-off parameter e decreases the accuracy of the approximation. 


Log[<5P] Log[<SP] 




(a) 


(b) 


FIG. 4: The relative error (SV) of the approximated to the numerically evaluated the 
probability, for k = 1 , e = 10~ 2 with (a) p = 0.1 and (b) p = 0.01 . 


Appendix B: Added-up probability 

The method of added-up probabilities, introduced by Audretsch and Spangehl in ref. |22j, 
gives a nice way of interpreting in-in amplitudes in terms of quantities, the definition of which 
is invariant to the nature of the out state. The problem to which it offers a solution is the 
following: in a setting with interacting fields on a dynamical background, when calculating 
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scattering amplitudes in the S-matrix approach, the effects of the interaction can not, in 
general, be separated from the cosmological pair creation. 

In a quasi-flat spacetime, where the cosmological particle production is negligible, one 
can do the usual Minkowskian perturbation theory. On the other hand, for strong grav¬ 
itational fields it is not clear whether and how one can neglect or subtract the effects of 
the cosmological particle production. An alternative possibility is to work with a notion of 
global measuring apparatus as in ref. [16], and interpret the in-in probabilities as the sole 
contribution of the mutual interaction, separate from the pure cosmological pair production. 

In the particular case when one of the fields is conformal, for example the Maxwell 
held in QED, the definition of the vacuum is unambiguous. This means that quanta of 
this held cannot be created freely from the vacuum, which in turn means that they are 
good indicators for the effects of the interaction, as the authors of ref. E2] point out. The 
physically measurable quantity then, they argue, is the probability of measuring a certain 
hnal state for the conformal held, regardless of the state of the other helds, given a hxed 
initial configuration. In QFT language this means that the summed probability is obtained 
by summing and integrating over the complete Fock space of the other helds in the hnal 
state. If we want for example, in the context of tree-level scalar QED, to calculate the 
summed probability for finding no photons, given a photon in the initial state, the summed 
probability is: 


'<4#(p.p'.*0 = £ I (Om; #|5 (1> |i ( s, a) ; op 


(Bi) 


# 


O1/3 1 *5'^ M 0m; #} ( Om; # | S 1 ' ) | 0^,), 

# 


where Om and 0^, represent the vacuum states of the Maxwell and scalar held, and # stands 
for the complete out Fock space for the scalar held. 


22 = f d 3 pi + f f d 3 pid 3 p 2 + - = 53 / " ’ / II d3pi 

# d J J n J J i = Q 


(B2) 


We observe by writing | 0 m,#) = | 0 m) 0 I#) i n (Bl), that we have obtained an identity 
operator for the scalar sector. We can thus insert any orthonormal base in place of #, in 
particular we can insert an in base. 


£i#x#l £ |m)(m| (B3) 

# in 
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This turns out to be the best solution for calculating such probabilities, as the summed 


probability turns into a finite sum of in-in probabilities. In the particular case of (Bl), the 
only term that contributes to the sum is the pair production process (fig. [I]). 

2 


W 


, add 
7 -># 


( l(p), l(p') I l(fc,A) ) 


(B4) 


This means that the in-in probability for one-photon pair production is equivalent, through 
the notion of the summed probability, with the probability that a photon is absorbed, 
regardless of the mechanism (regardless of ’’into what” it is absorbed). 

It may seem that this is not saying much, but we must note that the second part of the 
above statement describes a quantity which makes no assumption on the nature of the out 
state of the scalar field. 

As a further example, we evaluate the summed probability that is correlate to the one- 
photon annihilation process, the time-inversed version of fig. ([Tj) . This is interpreted as the 
probability of measuring a photon with momentum k in the out state (irrespective of the 
state of the scalar field), given a pair of scalar particles in the initial state. 

2 


add 


| ( 1 2 3 (fc,A)’ #I 5 ( 1 ) I°m; !(?), V)) 

# 

( l(jfc,A) I ^ 1 (p)> I(p') ) 

( l(p-fc)’ I(p') I 'I l(p)> l(p') ) 

( l(fc,A)> 1 (p)> l(p'-fc) I ^ ^1 -*■(?)> I(p') ) 

I d V ( l(fc,A)’ -*■(?")’ ^(- p"-k)i ■*■(?)> I(p') I ^ ^ W') ) 


(B5) 


The first term represents the one-photon annihilation process. The others are emission and 
triplet production processes, respectively, where the other particles pass through to the final 
state unchanged. These are the processes that are indistinguishable from the point of view 
of a photon counter that measures the number of photons in the final state [22]. 
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